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ABSTRACT 

Measurements  of  the  characteristic  time  spread  L  (in  seconds  )  of  surface  forward 
scattering  were  made  from  the  research  platform  FLIP  in  January  1992,  These  data 
were  compared  with  a  model  for  L  and  matched  the  model  to  within  an  average 
relative  error  of  30%.  The  inverse,  (in  hertz),  is  a  measure  of  the  coherence 
bandwidth  of  the  channel  and  therefore  important  in  both  designing  waveforms  and 
determining  signal-processing  strategies.  Using  the  model  for  L,  simulators  can  assess 
realistic  magnitudes  of  the  coherence  bandwidth  of  the  surface  channel  as  a  function 
of  source /receiver  geometry  and  environmental  conditions. 
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EXECUTIVE  SUMMARY 


The  Applied  Physics  Laboratory  at  the  University  of  Washington  (APL-UW) 
conducted  forward- scattering  measurements  between  20  and  50  kHz  from  the  re¬ 
search  platform  FLIP  oS  the  coast  of  California  in  January  1992.  This  work  was 
sponsored  by  the  Office  of  Naval  Technology  (now  ONR-T)  Torpedo  Environments 
Program  with  technical  management  by  the  Applied  Research  Laboratory,  Pennsyl¬ 
vania  State  University.  The  main  objective  wcis  to  obtain  coordinated  measurements 
of  acoustic  surface  scattering,  the  subsurface  bubble  field,  and  the  sea  state  in  order 
to  (1)  understand  more  fully  conditions  imposed  by  the  near-surface  environment 
that  affect  the  performance  of  torpedo  guidance  and  control  algorithms,  and  (2)  im¬ 
prove  high-frequency  acoustic  models  used  in  designing,  simulating,  and  predicting 
the  performance  of  underwater  weapon  systems. 

This  report  presents  estimates  of  the  characteristic  time  spread  L  (in  seconds) 
in  the  surface  forward  scattering.  The  inverse,  L~^  (in  hertz),  is  a  measure  of  the 
coherence  bandwidth  of  the  channel  and  therefore  important  in  both  designing  wave¬ 
forms  and  determining  signal-processing  strategies.  The  model  provides  a  means  to 
obtain  realistic  estimates  of  the  coherence  bandwidth  for  the  channel  as  a  function  of 
both  acquisition  geometry  and  environment. 

Fifty-eight  estimates  of  L  were  made  for  frequencies  between  20  and  50  kHz. 
ranges  of  540-1000  m,  grazing  angles  of  4.5°-17.6'’  and  wind  speeds  from  0-8  m/s 
(16  knots).  These  estimates  were  compared  with  a  model  for  L  and  matched  the 
model  to  within  an  average  relative  error  of  ±30%.  When  the  coherence  bandwidth 
L~^  is  estimated  with  the  model,  the  relative  error  translates  to  approximately  ±40% 
and  —20%. 

The  model  for  time  spread  is 


L  = 


2rir2  tan^o 
ri  ±  r2  c 


1 


where  ri  and  r2  are  the  incident  and  scattered  slant  ranges  (in  meters),  respectively, 
along  the  specular  path,  c  is  the  sound  speed  (in  meters  per  second),  and  0  is  the 
nominal  grazing  angle  of  the  surface  bounce  path.  The  model  requires  a  submodel 
for  the  rms  slope  s  =  tan  7oi  where 

s  =  ^0.0046  log,(2.lU2)  ,  U>lm 
s  =  0.0584  ,  C/  <  1  m  , 


and  U  is  the  wind  speed  (in  meters  per  second) 


measured  10  m  above  the  sea  surface. 


TR  9405  vi 


UNIVERSITY  OF  WASHINGTON*  APPLIED  PHYSICS  LABORATORY, 


The  coherence  bandwidth  estimates  as  determined  by  the  time- spread  estimates 
were  also  compared  with  a  limited  number  of  direct  estimates  of  the  frequency- 
coherence  function.  The  latter,  however,  were  encumbered  by  an  interference  between 
direct  and  surface-scattered  paths  and  by  an  interference  due  to  time  delays  between 
different  frequency  transmissions,  both  of  which  reduced  the  frequency-coherence 
function  at  a  frequency  separation  of  60-80  Hz. 
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1.  INTRODUCTION 

This  report  presents  results  of  an  experiment  to  measure  the  characteristic  time 
spread  L  (in  seconds)  in  surface  forward-scattered  intensity  at  frequencies  between  20 
and  50  kHz.  The  inverse,  L~^  (in  hertz),  is  a  measure  of  the  width  of  the  frequency- 
coherence  function  governing  frequency-selective  fading  and  is  therefore  important 
in  both  waveform  design  and  signal-processing  strategies.  The  estimates  of  L  are 
compared  to  a  model  for  L  based  on  geometric  variables  (range,  source  and  receiver 
depths,  and  surface  grazing  angle)  and  environmental  variables  (wind  speed  and  rms 
sea  slope).  The  model  provides  a  means  to  obtain  realistic  estimates  of  the  coherence 
bandwidth  for  the  channel  as  a  function  of  both  acquisition  geometry  and  environ¬ 
ment. 

The  measurements  were  made  from  the  research  platform  FLIP  off  the  coast  of 
California  in  January  1992.  The  objective  was  to  obtain  coordinated  measurements 
of  acoustic  surface  scattering,  the  subsurface  bubble  field,  and  the  sea  state.  The 
measurements  wdll  be  used  to  (1)  understand  more  fully  conditions  imposed  by  the 
near-surface  environment  that  affect  the  performance  of  torpedo  guidance  and  control 
algorithms,  and  (2)  improve  high-frequency  acoustic  models  used  in  the  modeling  and 
analysis  of  underwater  weapon  systems.  Additional  details  on  the  experiment  are 
given  in  Ref.  1,  which  also  summarizes  measurements  of  and  a  semiempirical  model 
for  the  mean  energy  loss  in  surface-bounce  paths  due  to  extinction  from  near-surface 
bubbles. 

Preliminary  remarks  and  a  summary  of  notation  used  in  this  report  are  given 
in  Section  2.  The  experimental  geometry  and  results  of  the  time-spread  analysis 
and  environmental  sea-state  measurements  are  given  in  Section  3,  along  with  model 
comparisons.  Five  signal  suites  were  also  amenable  to  estimation  of  the  actual  shape 
of  the  frequency-coherence  function;  although  this  data  set  is  sparse,  the  results  are 
included  for  comparison  in  Section  4.  A  discussion  and  summary  are  presented  in 
Section  5. 
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2.  PRELIMINARY  REMARKS  AND  NOTATION 

A  brief  summary  of  notation  and  usage  is  presented  here.  Further  details  on  this 
subject  are  available  in  Refs.  2-5  and  in  the  additional  references  cited  below. 

The  scattering  function  i2j(T,  4>)  is  a  measure  of  the  average  delay,  or  time  spread. 
(t)  and  frequency  spread  {(f>)  that  a  transmitted  signal  will  be  subjected  to  by  the 
channel  (on  an  individual  ping  basis,  this  function  is  known  as  the  spreading  function 
and  varies  from  ping  to  ping).  For  the  forward-scattering  measurements  made  from 
FLIP,  the  channel  is  defined  as  a  single  interaction  with  a  time-varying  random  sea 
surface.  Such  a  channel  is  usually  time  dispersive,  because  there  are  multipath  arrivals 
leading  to  a  spread  of  delays,  and  frequency  dispersive,  because  the  medium  is  time- 
varying,  causing  a  broadening  of  the  received  signal  spectrum  (or  Doppler  spread). 
Both  time  and  frequency  dispersion  render  the  channel  doubly  spread®;  the  property 
of  angular  dispersion  adds  yet  another  dimension  of  channel  spreading. 

In  a  manner  analogous  to  the  joint  probability  density  function  for  two  random 
variables,  in  this  case  r  and  (f),  one  can  define  a  marginal- scattering  function  for 
which  one  of  the  variables  has  been  integrated  out.®  We  define  L  (in  seconds)  as 
the  characteristic  width  for  a  marginal-scattering  function  linked  to  time  spread  and 
B  (in  hertz)  as  the  characteristic  width  for  a  marginal-scattering  function  linked 
to  frequency  spread.  More  precisely,  we  can  define  L  and  B  .ts  the  square  root  of 
the  second  moment  of  their  respective  marginal-scattering  functions.  Under  some 
conditions,  the  scattering  function  may  be  separable'  such  that 

R,{t,4>)  =  Rr{T)RA0)  ^  (1) 

where  Rt{t)  is  the  so-called  delay,  or  time-spread,  profile  and  Roio)  is  the  Doppler, 
or  frequency-spread,  profile. 

If  scattering  from  the  sea  surface  is  incoherent,  then  the  intensity  time-spread 
function  and  frequency-coherence  function  are  related  by  Fourier  transforms.®  Ac¬ 
cordingly,  the  inverse,  L~^  (in  hertz),  is  an  approximate  measure  of  the  coherence 
Viandwidth  of  the  channel.  In  other  words,  transmitting  signals  with  a  bandwidth 
W  >  L~^  may  result  in  frequency- selective  fading,  while  for  W'  <  L~^  all  frequencies 
behave  more  or  less  the  same  in  terms  of  fading  statistics. 

For  completeness,  an  estimate  for  B  can  be  derived  using^ 

B  =  0.mUfsme  ,  (2) 

where  /  is  the  frequency  in  kilohertz,  U  is  the  wind  speed  in  meters  per  second 
measured  10  m  above  the  sea  surface,  and  6  is  the  sea-surface  grazing  angle.  Based 
on  Eq.  (2),  B  is  approximately  5  Hz  for  the  range  of  grazing  angles,  wind  speeds,  and 
frequencies  that  apply  to  the  FLIP  forward-scattering  measurements. 
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3.  EXPERIMENTAL  GEOMETRY,  TIME-SPREAD  MEASURE¬ 
MENTS,  SEA-SLOPE  COMPUTATIONS,  AND  MODEL 
COMPARISONS 

The  geometry  used  to  measure  the  time  spread  in  surface  forward-scattered 
intensity  is  shown  in  Figure  1.  Signals  were  transmitted  from  one  of  three  ITC-1032 
transducers  suspended  at  depths  of  27,  57,  and  147  m  from  the  spar  buoy.  The  spar 
buoy  was  tethered  to  FLIP  by  a  1000-m-long  cable  and  ranged  from  500  to  1000  m 
from  the  vessel.  The  receiver  options  at  FL/P  consisted  of  the  lTC-1032  transducer 
located  on  the  tip  of  the  12-m-iong  subsurface  boom  attached  to  FLIP'S  hull  at  a 
depth  of  28.5  m,  the  line  array,  and  the  Mk-46  array.  Both  arrays  were  attached  to 
flip's  hull  at  a  depth  of  66  m. 

Frequency  spectra  for  sea-surface  displacement  Sr,,j{f).  made  with  a  vertical- 
incidence  scatterometer,  were  taken  every  20  minutes  by  A.  Jessup  of  APL-l'W  . 
From  the  displacement  spectra,  we  computed  a  large-scale  mean-square  surface  slope 
using 


Figure  1.  Experimental  geometry  for  measurements  of  time  spread  and  frequency 
coherence  in  surface  forward  scattering.  Measurements  were  made  using 
all  three  sources  suspended  from  the  spar  buoy  and  all  three  receiver  systems 
attached  to  FLIPs  hull.  The  spar's  buoy  range  varied  between  500  and 
1000  m. 
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4' =  ^  >  (3) 

5  Jo 

where  p  is  a  gravitational  constant  and  fc  is  the  integration  cutoff  frequency.  For  a  sea- 
surface  wavenumber  spectrum,  numerical  simulations  by  McDanieF“  suggest  setting 
a  wavenumber  cutoff  equal  to  0.4  times  the  magnitude  of  the  acoustic  wavenumber 
k.  Upon  using  the  deep-water  dispersion  relation,  this  gives  fc  equal  to  sj if) Agk) ! '1~ . 
The  spot  size  of  the  scatterometer  prevents  frequencies  ^  0.6  Hz  from  being  ade¬ 
quately  resolved.  To  correct  for  this,  an  /~®  tail  is  appended  to  S„„  for  frequencies 
>  0.6  Hz."  The  rms  slope  shs  significant  wave  height  (SWH),  and  wind  speed  for  each 
acoustic  measurement  are  listed  in  Table  1.  Additional  details  on  the  environmental 
measurements  are  discussed  in  Refs.  1  and  12. 

We  estimate  L  in  the  following  manner.  A  modeP^  for  the  intensity  impulse 
response  is  numerically  convolved  with  the  particular  transmit  pulse  envelope  used 
in  the  FLIP  measurements  (t.3'pically  1-8  ms  cw').  The  result  gives  a  model  for  the 
time  dependence  of  surface-scattered  intensity  which  includes  /,  as  a  free  parameter. 
Wo  then  estimate  L  by  matching  this  model  with  ensemble- averaged  data  in  a  least 
squares  sense,  such  as  in  the  example  shown  in  Figure  2  where  L  =  3  ms.  Results  for 
58  such  estimates  made  for  frequencies  between  20  and  50  kHz.  ranges  of  500-1000  m. 
grazing  angles  of  and  wind  speeds  of  0-8  m/s  are  listed  in  Table  1.  In  the 

table,  SD  =  source  depth,  RD  =  receiver  depth,  and  T  =  pulse  length.  The  variables 
sa  and  Tmod  are  discussed  in  the  next  section. 


Figure  2.  Model  for  the  time  dependence  of  surface-scattered  intensity  (thin  line) 
fitted  to  data  (circles)  using  L  =  3  ms. 
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Table  1.  Summary  of  time-spread  database  derived  from  FLIP  experiment. 
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4.  MODEL  VS  DATA  COMPARISONS 

A  model  for  the  characteristic  time  spread  L  in  forward  scatter 

r,  +  ra  c 

where  ri  and  ra  are  the  incident  and  scattered  slant  ranges  (in  meters),  respectively, 
along  the  specular  path,  c  is  the  nominal  sound  speed  {in  meters  per  second),  and  6  is 
the  nominal  grazing  angle  of  the  surface  bounce  path.  1  model  requires  a  submodel 
for  the  rms  slope  s  =  tan  70.  For  this,  we  use  a  model  baised  on  a  logarithmic  fit*^  of 
the  classic  set  of  optical  glitter  measurements  by  Cox  and  Munk"  of  the  mean  square 
slope  for  a  slick-covered  surface, 

s  =  Y^0.0046log,(2.1tf2)  ,  U>\m 
s  =  0.0584  ,  <  1  m  .  (5) 

The  logarithmic  wind  speed  dependence  in  Eq.  (5)  from  Ref.  15  is  restricted  to  wind 
speeds  U  >  I  m/s.  For  wind  speeds  less  than  1  m/s,  we  set  s  equal  to  0.0584 
(equivalent  to  Eq.  (5)  evaluated  at  1  m/s)  because  data  on  sea  surface  slope  for 
extremely  low'  wind  speeds  are  too  scattered  to  specify  a  more  precise  functional 
form. 

The  time  spreads  Tmod  predicted  using  the  model  in  Eq.  (4)  and  the  submodel  in 
Eq.  (5)  are  given  alongside  each  measured  time  spread  L  in  Table  1.  The  predictions 
have  a  relative  uncertainty  of  approximately  ±30%.  The  rms  value  of  the  model 
(Tmod)  minus  data  (X)  vectors  is  0.75  ms.  We  conclude  that  Eq.  (4),  along  with  the 
submodel  for  rms  sea  slope,  enables  one  to  assess  realistic  magnitudes  of  the  coherence 
bandwidth  (defined  as  L~^)  as  a  function  of  both  acquisition  geometry  (ri,  r2,  and 
6)  and  environment  (U).  The  relative  uncertainty  in  L  translates  to  approximately 
±40%  and  —20%  for  coherence  bandwidth  predictions. 

We  can  also  infer  an  estimate  of  the  large-scale  sea  slope  s  from  the  original 
measurements  of  L  and  the  geometric  variables  listed  in  Table  1.  To  do  this,  Eq.  (4) 
is  solved  numerically  using  the  Newton- Raphson  method’^  to  find  the  unknown  lo- 
The  sea-slope  estimates  (equal  to  tan  70)  derived  from  the  acoustic  estimates  of  time 
spread  are  referred  to  as  in  Table  1  and  are  plotted  vs  wind  speed  in  Figure  3 
along  with  sw  and  s  computed  using  Eq.  (5).  The  scatter  show’n  in  Figure  3  is  quite 
typical  of  slope- vs- wind-speed  data.  The  model  for  s  based  on  the  Cox  and  Munk 

‘The  slope  measurements  by  Cox  and  Munk  were  made  in  both  clean  water  surface  conditions 
an  il-slick-covered  surface  conditions.  The  effect  of  the  oil  slick  was  to  effectively  damp  out  surface 
waves  with  a  wavelength  less  than  about  30  cm.  Therefore  we  use  the  slick-covered  measurements 
because  they  are  more  representative  of  slope  conditions  sampled  by  finite  acoustic  wavelengths 
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Figure  3.  Model  for  rms  large-scale  sea  slope,  Eq.  (5),  vs  slope  estimates  derived  from 
the  time-spread  database  and  inversion  of  Eq.  (4)  (+)  and  slope  estimates 
derived  from  wave-gauge  data  and  Eq.  (3)  (o). 


data,  Eq.(5),  describes  the  sa  data  reasonably  well  and  tends  to  form  an  upper  bound 
to  the  sw  data.  Recall  that  the  sw  estimates  depend  on  an  integration  cutoff  that 
increases  with  acoustic  frequency.  Thus  the  multifrequency  runs  (20,  30,  40,  and 
50  kHz)  listed  in  Table  1  will  show  up  as  a  clustering  of  points  in  Figure  3,  with  the 
highest  frequency  giving  a  shghtly  higher  estimate  of  sw.  Furthermore  the  f~^  tail, 
appended  to  Sr,T,  for  frequencies  >0.6  Hz,  influences  the  integral  in  Eq.  (3).  A  tail  of 
for  example,  would  increeise  the  final  estimates  for  as  derived  from  Eq.  (3). 


TR9405  7 


UNIVERSITY  OF  WASHINGTON*  APPLIED  PHYSICS  LABORATORY. 


5.  FREQUENCY  COHERENCE  FUNCTION  ESTIMATES 

Estimates  of  the  frequency-coherence  function  for  forward  scattering  were  made 
using  signal  suites  specified  by  Dr.  D.  Ricker  of  ARL-PSU;  these  suites  are  referred  to 
as  Arctic,  Costas-1,  and  Costeis-2.  The  received  data  from  the  Arctic  pulse  are  stored 
in  a  file  called  FWDmk46.101.  The  original  (transmitted)  pulse  is  560  ms  long  and 
constitutes  seven  frequency  chips,  with  each  chip  having  a  different  frequency  offset 
from  the  carrier  frequency.  The  received  data  from  the  Costas-1  pulse  are  stored 
in  the  file  FWDmk46.102.  The  original  pulse  is  500  ms  long  and  divided  into  three 
periods;  the  first  200  ms  contains  27  frequency  chips,  the  following  100  ms  is  silence, 
and  the  last  200  ms  contains  27  frequency  chips.  The  received  data  from  the  Costas-2 
pulse  are  stored  in  the  files  FWDmk46.103  and  FWDmk46.105.  Each  of  the  above 
files  consists  of  25  pings  recorded  on  four  channels.  For  the  Costas-1  and  Costas-2 
pulses,  the  received  signal  is  split  into  two  parts  (see  Table  2)  because  of  the  100-ms 
silence  between  the  two  transmitted  waveforms:  each  part  is  the  channel  response  to 
waveform  A,  B,  or  C.  The  offset  frequencies  for  the  various  chips  used  in  waveform  A 
are  listed  in  Table  3.  Waveform  B  is  the  same  as  A  with  the  sign  reversed;  waveform 
C  is  a  200-ms  cw  pulse  with  zero  offset  from  the  carrier  frequency. 
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Frequency  (Hz) 

-498.1 

-273.1 

-203.8 

-411.5 

-376.9 

-290.4 

-134.6 

-359.6 

-550 

-446.2 

-117.3 

-186.5 

-325 

-307.7 

-342.3 

-221.2 

-169.2 

•515.4 

-151.9 

-463.5 

-255.8 

-428.8 

-532.7 

-100 

-394.2 

-480.8 

-238.5 

Table  3. 

Offset  frequency  for  the  27  chips  in  Costas 
waveform  A.  The  actual  transmitted  fre¬ 
quency  equals  a  carrier  frequency  plus  that 
shown  in  the  table.  The  order  of  chip  trans¬ 
mission  is  shown  from  top  to  bottom;  e.g., 
the  first  entry  f— 498.1  Hz)  was  the  first 
chip  to  be  transmitted. 


The  received  signal  is  processed  through  a  filter  bank,  where  the  transfer  function 
of  each  cell  is  weighted  with  a  Hanning  window  with  a  central  frequency  and  width 
equal  to  those  of  the  corresponding  frequency  chip.  Each  filter-bank  output  is  cross 
correlated  with  every  other  output.  The  coherence  between  the  responses  of  two 
frequency  chips  is  defined  as 


(y/.(0y/,(0) 


(6) 


where  /,  and  fj  are  the  central  frequencies  of  the  and  the  chips,  and  yj,{t)  is 
the  output  of  the  cell  of  the  filter  bank.  The  coherence  magnitude  corresponds 
to  |7ij(T  =  0)1  and  indicates  the  amount  of  coherence  between  the  chip  centered  at 
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frequency  /,  and  that  centered  at  frequency  fj.  The  final  result  is  an  N  x  N  matrix 
of  l7ij(^  =  0)1'  where  N  is  the  number  of  chips;  an  example  is  shown  in  Figure  ■!. 

The  chips  can  be  ordered,  or  numbered,  by  their  frequency  offset  or  by  their  time 
of  transmission  within  the  waveform.  In  Figure  4  the  bottom  and  left  axes  show  the 
chip  numbers  according  to  offset  frequency  sequencing,  and  the  top  and  right  axes 
show  the  chip  numbers  according  to  the  transmission  time  sequencing.  For  example, 
the  chip  with  an  offset  frequency  of  -151.9  Hz,  has  a  frequency  sequence  number  of  24 
and  a  time  sequence  number  of  19  (Table  3).  The  best  coherence  is  achievable  when 
chips  are  close  in  both  frequency  and  time,  e.g.,  two  chips  that  are  close  in  frequency 
but  far  apart  in  time  will  still  show  weak  coherence.  For  example,  in  Figure  4  the 
neighboring  frequency  chips  4  (-498.1  Hz)  and  5  (-480.8  Hz)  show  weak  coherence 
because  they  were  transmitted  far  apart  in  time. 

Figure  4  summarizes  all  the  coherence  estimates  for  one  run,  but  it  is  admittedly 
difficult  to  interpret.  A  simpler  view  is  shown  in  Figures  5-9  which  show  the  coherence 
of  pairs  of  temporally  neighboring  chips  as  a  function  of  their  frequency  separation. 
Note  that  the  analysis  was  performed  using  the  ensemble  average  (25  pings)  for  one 
quad  and  for  the  coherent  sum  of  four  quads  of  the  Mk  46  array,  and  no  significant 
difference  was  observed  in  the  results. 
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Figure  5.  Magnitude  of  frequency  coherence  vs  frequency  difference  A  /  based  on  tem¬ 
porally  neighboring  chips  for  file  FWDmk46.101  and  pulse  type  Arctic. 


Figure  6.  Magnitude  of  frequency  coherence  vs  frequency  difference  A/  based  on  tem¬ 
porally  neighboring  chips  for  file  FWDmk46.102  and  pulse  type  Costas- 1  A. 
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Figure  7.  Magnitude  of  frequency  coherence  vs  frequency  difference  Af  based  on  tem¬ 
porally  neighboring  chips  for  file  FWDmk46.102  and  pulse  type  Costas- 1  B. 


Figure  8.  Magnitude  of  frequency  coherence  vs  frequency  difference  Af  based  on  tem¬ 
porally  neighboring  chips  for  file  FWDmkfS.  103  and  pulse  type  Costas-2  A. 


Figure  9.  Magnitude  of  frequency  coherence  vs  frequency  difference  Af  based  on  tem¬ 
porally  neighboring  chips  for  file  FWDmk46.I05  and  pulse  type  Costas-2  A. 
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6.  DISCUSSION  AND  SUMMARY 

The  acoustic  runs  for  which  the  frequency  coherence  was  estimated  and  shown 
in  Figures  5-9  were  taken  within  1  hour  of  each  other.  The  source  and  receiver  depths 
were  57  and  66  m,  respectively;  the  range  was  676  m,  the  grazing  angle  10°,  and  the 
wind  speed  4  m/s.  Using  these  variables  in  Eqs.  (4)  and  (5),  we  would  expect  the 
time  spread  L  to  be  about  2.5  ms  and  therefore  the  coherence  bandwidth  to  be  about 
400  Hz.  The  coherence  bandwidth  estimated  from  Figures  5-9  is  clearly  less  than 
400  Hz,  more  on  the  order  of  50  Hz  (full  bandwidth). 

A  probable  reason  for  this  difference  is  that,  for  these  runs,  we  were  unable 
to  effectively  separate  the  direct  path  from  the  surface-reflected  path  which  arrived 
~7  ms  later.  Thus  the  frequency-coherence  functions  shown  in  Figures  5-9  suffer  from 
interference  between  the  direct  and  surface-bounce  paths  (although  the  frequency 
spacing  in  Figure  5  is  not  small  enough  to  show  this).  Note  also  that  there  is  a  time 
delay  of  7.5  ms  between  neighboring  frequency  chips  within  the  Costas  waveforms  used 
for  the  coherence  estimates  in  Figures  6-9.  The  combination  of  direct  and  surface 
interference  and  the  time  delay  between  neighboring  frequency  chips  significantly 
reduces  the  coherence  for  frequency  differences  A/  satisfying  relation  |A/|AT  ~ 
0.5,  where  AT  i?  either  the  travel-time  difference  or  the  chip  time  delay.  Frequency 
differences  j  A  f\  in  the  neighborhood  of  60  to  80  Hz  meet  this  criterion,  with  this  effect 
shown  in  an  approximate  sense  in  Figures  6~9  (though  less  clearly  in  Figure  7).  Note, 
however,  that  in  these  same  figures  the  coherence  rises  for  |A/j  sufficiently  greater 
than  about  80  Hz  but  less  than  about  200  Hz.  This  is  reasonably  consistent  with  a 
coherence  full  bandwidth  of  ~400  Hz,  as  predicted  by  the  time-spread  estimates. 

In  summary,  our  main  findings  pertain  to  the  parameter  L,  defined  as  the  char¬ 
acteristic  time  spread,  in  seconds,  in  the  surface  forward- scattered  intensity.  The  58 
estimates  of  L  were  compared  with  predicted  values  based  on  a  model,  Eqs.  (4)  and 
(5),  and  the  model  matched  the  data  within  a  relative  uncertainty  of  ±30%.  The 
inverse,  L~^  (in  hertz),  is  an  approximate  measure  of  the  coherence  bandwidth  of 
the  channel.  Using  the  model  for  L,  simulators  can  assess  realistic  magnitudes  of  the 
surface-channel  coherence  bandwidth  as  a  function  of  source /receiver  geometry  and 
environmental  conditions. 

We  were  unable  to  match  estimates  of  the  frequency  coherence  function  with 
the  time-spread  data  because  of  interference  between  direct  and  surface-bounce  paths 
and  interference  due  to  time  delays  between  different  frequency  transmissions,  both 
of  which  reduced  the  magnitude  of  the  frequency-coherence  function  for  |A/|  in  the 
range  60-80  Hz.  Coherence  properties  away  from  this  range  were  reasonably  consis¬ 
tent  with  the  expected  coherence  bandwidth  based  on  the  time-spread  meeisurements. 
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